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Flat magnetic nano-elements are an essential component of current and future spintronic devices. 

By shaping an element it is possible to select and stabilize chosen metastable magnetic states, 
control its magnetization dynamics. Here, using a recent significant development in mathematics of 
conformal mapping, complex variable based approach to the description of magnetic states in planar 
nano-elements is extended to the case when elements are multiply-connected (that is, contain holes 
or magnetic anti-dots). We show that presence of holes implies a certain restriction on the set of 
magnetic states of nano-element. 

PACS numbers: 75.60.Ch, 75.70.Kw, 85.70.Kh 


While existence of topological solitonsi as metastable 
states in infinite 2-d ferromagnets was known theoreti¬ 
cally for quite a long timej2r— real systems have finite 
size. Their boundary imposes restrictions on the set of 
low-energy topological states, making it equivalent to the 
set of rational functions of complex variable with real 
coefficients^ as opposed to complex coefficients in the 
laterally unconstrained thin film case: 2 Also, restricted 
geometry implies the possibility of formation of half¬ 
vortex states, pinned at the side of the planar magnet^ 
which are topologically similar to the boundary states 
in the fractional quantum Hall systems and topological 
insulators. The shape of nano-element enters this com¬ 
plex description^ via its conformal mapping to the unit 
disk. It is also possible to express Landau-Lifshitz mag¬ 
netization dynamics directly in terms of these complex 
functions using Lagrangian approach. 7 ' The complex de¬ 
scription of magnetization states^ and their dynamics^ 
constitutes a complete set of analytical tools to study 
magnetic textures in simply-connected nanomagnets far 
from magnetic saturation. Here this description is ex¬ 
tended for the case when planar nano-magnet is multiply- 
connected. 

Let us briefly introduce the description for simply- 
connected case£ In small enough magnets the surface ef¬ 
fects dominate the volume ones, and also the exchange in¬ 
teraction is more important than the magnetostatic one. 
The former is typical for any small systems, while the 
latter is easy to understand based on the representation 
of the magnetostatic energy in terms of the interaction 
of fictitious magnetic charges of opposite signs. The to¬ 
tal magnetic charge is always zero, thus, whatever the 
distribution of the magnetization is, as the size of the 
magnet decreases, the positive and negative charges are 
brought closer together, so that their positive self-energy 
gets more and more compensated by their negative mu¬ 
tual interaction energy. The scaling of the exchange en¬ 
ergy of the set of such an arbitrary magnetization distri¬ 
butions follows the volume of the magnet and does not 
have such an additional reduction. Thus, in small mag¬ 


nets the exchange interaction becomes more important. 
To make use of this energy hierarchy, we build an approx¬ 
imate expressions for magnetic states by minimizing the 
energy terms sequentially (as opposed to minimization 
of their sum, which would result in the exact theory). 
The process of sequential minimization is analogous to 
sieving. 

First, consider all possible magnetization distributions 
with fixed length of magnetization vector \m\ = 1: 


mx + imy = 


m z = 


2 w(z,z) 

1 + w(z,~z)w(z, zc) 
1 — w(z,z)w(z, zc) 

1 — w(z,z)w(z, zc) ’ 


(1) 


where fh = {mx , my, mz} = M/Ms is the local magne¬ 
tization vector, expressed in units of saturation magne¬ 
tization Mg ; z = X + iY , with X, Y and Z being the 
Cartesian coordinates (the element is assumed to be the 
a flat cylinder with axis, parallel to Z and magnetization 
distributions are assumed to be ^-independent), w(z,z) 
is a function (not necessarily meromorphic) of complex 
variable 2 . 

Second, out of all these functions select the ones, 
which, additionally to having the constant length of the 
magnetization, minimize the exchange energy (these are 
the famous Belavin and Polyakov solitons^). Then fur¬ 
ther restrict the set of the remaining functions twice by 
selecting those, which give minimum for the energies of 
face magnetic charges and those, which minimize (in fact, 
totally avoid) the side magnetic charges. The final result^ 
of such a selection will be the following representation for 
w(z,z) 


f f(z)/e l_ \f{z)\ < ei 

w(z,z) = < f{z)/\Jf{z)J(z) ei< |/0)| <e 2 , (2) 

{ f(z)/e 2 |/0)| >e 2 

where ei and e 2 are real positive constants. The complex 
function f(z) is a solution of Riemann-Hilbert boundary 
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FIG. 1. Equilibrium and transient magnetization states in 
ferromagnetic nano-disk following from the Eq. [3] a) centered 
magnetic vortex; b) “leaf” state; c) displaced magnetic vortex; 
d) “C”-like state. 


value problem of finding the meromorphic function in 
the domain V , which corresponds to the face of the pla¬ 
nar nanoelement, having no normal components to the 
domain’s boundary. This problem usually has many so¬ 
lutions. For simply-connected case their set is equivalent 
to the set of rational functions with real coefficients (con¬ 
stants), whose zeros (poles) correspond to vortex(anti- 
vortex) centers £ 

For example, in a unit disk the subset of states with 
no antivortices can be expressed^ as 

f disk (z)=izc + A-Az 2 , (3) 

where c and A are a real and a complex constants respec¬ 
tively. While the equilibrium values of these constants 
can be found only by minimizing the total (including the 
energy of the volume magnetic charges) energy of the 
nanomagnet,— the expression J3]) can already be useful to 
pinpoint that there are several types of magnetic states in 
the disk, such as centered magnetic vortex (when A = 0), 
the so-called “leaf” state (when c = 0) and C-like magne¬ 
tization state (when 4AA > c). These states are shown 
in Figure 1. The solution ([3]) can be extended when there 
are additional vortex-antivortex pairs present^ allowing 
to describe both transient dynamical states as well as 
finite (simple and topologically charged domain walls). 

Before proceeding to the consideration of multiply- 
connected case, let us note a general property of the 
above Riemann-Hilbert boundary value problem. If it 
has a solution f(z) in the region z £ T>, the correspond¬ 
ing solution u(t) in another region t £ D', connected to 
T> via the conformal mapping t = T(z ), can be expressed 



FIG. 2. An example of a quadruply-connected circular do¬ 
main with an outer circular boundary Go and inner bound¬ 
aries Cj, j — 1, 2, 3. 


as 


= f(z)T'(z) | z=T( -i) (t ) , (4) 

where T^ _1 l(t) is the inverse of the conformal mapping 
t = T(z), which is always defined. This expression ap¬ 
plies both to simply- and multiply-connected cases and 
allows to express u(t) in an arbitrarily shaped cylinder, 
based on its expressions f(z) for some canonical domains. 
In the simply-connected case these canonical domains 
can be chosen, based on convenience only, to be the unit 
diskf£ the half-plane^ or any other simply connected re¬ 
gion for which the solution can be written explicitly. Se¬ 
lection of the canonical domain is trickier in the multiply- 
connected case, since not every pair of regions with the 
same connectivity can be conformally mapped into each 
other. 

Thanks to the Koebe’s theorems, it is well known that 
it is possible to define parametrized canonical families of 
multiply-connected regions, which, after parameter ad¬ 
justment, can be mapped to an arbitrary region of the 
same connectivity. Among them is the familj*^ of circu¬ 
lar domains with cut out inner circles shown in Fig. [2] 
Positions and radii of the circular holes in these domains 
are not arbitrary and are dictated by the shape of the 
target multiply-connected region V to which f(z) may 
be mapped. Because such canonical regions by them¬ 
selves have physical significance (such as the simplest 
case of concentric ring, which is doubly-connected), let 
us choose them as the basis and set the goal of express¬ 
ing the solutions f(z) of the Riemann-Hilbert problem in 
these regions. 

But first we need to address the question of exis¬ 
tence of such solutions and their general properties. For¬ 
mally, given a finitely connected flat domain T> with good 
enough boundary (say, piecewise analytic), the problem 
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is to find a meromorphic function f(z) in the domain, 
satisfying the boundary condition 

Re(f(z)n{z)) =0, z £ dV , (5) 

where n{z) is a normal to the boundary of the domain. 
The solution f(z) may have zeros and polynomial singu¬ 
larities at the boundary. This (Riemann-Hilbert) bound¬ 
ary condition can be reformulated as follows: 


of T> to V if the latter is a multiply-connected polyg¬ 
onal domain.— It also has a number of applications in 
problems of optimization and computation . 16,17 

Informally (for the formal definition see e.g. Ref. [bfh . 
the Schottky-Klein prime function w(z, C) for a specific 
multiply-connected circular domain V can be thought as 
a generalization of the difference 

z-C = wi(z,C)> (?) 


Im. 



z € dV, 


( 6 ) 


which in multiply connected case becomes 
w = w(z, 


( 8 ) 


which means that meromorphic differential (d z)/f(z) is 
real (i.e. the restriction of this differential to the bound¬ 
ary is a real differential). The latter condition is confor¬ 
mally invariant, which, in particular, implies (|d| . 

To describe the set of all real differentials in T >, the 
useful notion is the double of the domain. The latter is a 
compact Riemann surface X(T>) of genus g equal to the 
number of the boundary components of D minus one. 
It is made of two copies of the domain glued along its 
boundaries, z being the conformal coordinate on one copy 
and z — on the other. This surface admits natural an- 
ticonformal involution J (or reflection ), the interchange 
of copies, which is stationary exactly on the boundary 
components of the domain V. The reflection acts on the 
differentials and the condition for the differential di 7 to 
be real is exactly the following: 


In a simply-connected case products of such differences 
can be used to build the rational functions of complex 
variable, such as those equivalent to the set of Belavin- 
Polyakov solitonsS, or a similar set of states of finite nano- 
magnet^ or those, entering the Schwarz-Christoffel for¬ 
mula for the conformal mapping of polygons. General¬ 
ization of the latter to the multiply-connected case was 
done by Crowdy,— here we shall outline a generalization 
of the former. 

In the simplest doubly-connected case of a concentric 
ring with an outer radius of 1 and an inner radius of q < 1 
the Schottky-Klein prime function can be expressed as a 
product 


Wr{z,C) = {z-Q 


nr = i(^g 2 vo(i-g 2fc cA-) 

( nr = i ( i -<? 2fe )) 2 


(9) 


J dg = dg. 


or written via the q-Pochhammer symbols (a; q) n as 


Corollary: The following topological phenomenon 
can be observed: 


w r (z,C) 


z( (C/z,g 2 )oo (z/(,q 2 ) 00 

((A, A) oo ) 2 


( 10 ) 


tt{/(~) poles} — (({/A) zeroes} = boundaries} — 2 

where the poles and the zeros are counted with their mul¬ 
tiplicities and for those at the boundary the multiplicity 
should be divided by two. 

Proof: The solution f(z) corresponds to the mero¬ 
morphic differential dr? := (d z)/f(z) on the double of 
V. The degree of the divisor of a differential is 2g — 2. 

This implies that we can’t control independently num¬ 
ber of zeros and poles of f{z). 

The natural recipe to cook a real differential is sym- 
metrization: take any meromorphic differential dg on 
the surface, then its symmetrization dg + J dg is mero¬ 
morphic and real. This recipe is easy to use once we 
have a representation of the double X{V) as an algebraic 
curve. This representation, for the case of circular mul¬ 
tiply connected domains can be represented via a series 
over the elements of the corresponding Schottky group.— 
The evaluation of such series, however, can be very incon¬ 
venient if approached directly. That’s why in the follow¬ 
ing we will give the representation of the solution f(z) via 
the Schottky-Klein prime function , 12 ’ 13 which not only 
admits an efficient numerical evaluation,— but can also 
be directly used in building the conformal map t = T(z) 


This representation will be used in some of the following 
examples. 

To build the solutions of the Riemann-Hilbert problem 
with the specified positions of vortices (or anti-vortices) 
let us define two auxiliary functions, following from those, 
introduced in the Section 5 of Ref.[ 15f 


Fi{z, Ci> C2) 


AA.C) 


w(Z)Cl) 

w(z, C2)’ 
w(zX)w(z 7 1 /C) 
w{z,Ow{z,i/c)' 


( 11 ) 

( 12 ) 


which, provided Ci and C2 belong to the same inner 
or outer circle Cj (with j = 0 , 1 , 2 ,...) of the multiply- 
connected circular domain and C is any point inside of 
it, have constant argument (complex phase) on all the 
boundaries Cj of the said domain. This means that any 
product of these functions will have the constant argu¬ 
ment on the boundaries of the circular domain too, which 
implies that the logarithmic derivative of this product 


g(z) = vpA°g ( 


( 13 ) 










FIG. 3. Finite domain walls in a ring of q = 0.3, 
given by m and m with ei = 0 , e 2 —» oo, g(z) = 
(d/dz)\o S F 1 (z,e- 8 ™ /9 ,e™ /9 ) and f(z)=g(l/z). 



FIG. 4. A domain structure with several finite domain 
walls in a ring of q = 0.7, with g(z) containing now 
the product of three Fi functions: Fi(z, e 8,r ^ 9 , e _7rl ^ 9 ), 
F\(z, qe 4m / 3 , qe~ 2m/3 ) and Fi(z, qe™ /3 , qe~ nl/3 ). The rest of 
parameters are the same as in Figure [3] 


FIG. 5. Vortex domain walls in a ring of q = 0.6, 
with g(z) containing the product of Fi(z, e 4 ™/ 3 , e nl ^ 3 ) and 
F\(z,qe llnl ^ 9 ,qe 2vx ^ 9 ). The rest of parameters are the same 
as in Figure [3] 



FIG. 6 . Antivortex domain walls in a ring. The parameters 
for this plot are the same as in Figure [SJ except that now 
f(z) = g(z). 


will be an analytic function of 0 satisfying the condi¬ 
tion (0 at the boundary of the multiply-connected cir¬ 
cular domain. That is, f(z) = g{z) will be a solution 
of the Riemann-Hilbert problem with the specified (by 
the parameters Ci,m, C2 ,mi Cn) positions of zeros and 
poles. It will have some additional zeros and poles too, 
so that the previously mentioned topological constraint is 


always satisfied. Also, it is easy to show that the function 
f(z ) = < 7 ( 1 / 2 :) corresponds to a “mirror” solution of the 
problem, where the vortices, specified by the parameters 
£, correspond to the antivortices and vice versa. This 
is illustrated by several examples in Figures [3J7] some 
of which closely resemble well known magnetization tex¬ 
tures, observed in ferromagnetic nanorings. 

From these, the corresponding magnetization textures 
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FIG. 7. A vortex in a triply connected domain, with the cut 
out circles of the radii 0.2 and 0.3, located at x=0.5 and x=- 
0.5 respectively. The function g(z) = ( d/dz ) logF 2 (z, —0.52 j) 
and f(z) = g(z). The Schottky-Klein prime function was 
evaluated numerically, using the Matlab package developed 
in Ref. [l4|. 


in the nano-magnets of the same high connectivity but 
other shapes can be derived according to the Eq. [4] us¬ 
ing the conformal mapping (e.g. the mapping to the 
multiply-connected polygonal domains, which are also 
expressed in terms of Schottky-Klein prime function^). 

Thus, we have built an approximate analytical repre¬ 
sentation of the low-energy magnetization states in pla¬ 
nar multiply-connected nanomagnets. It is parametrized 
via the positions of some of vortices or antivortices in the 
magnet. Not all topological singularities can be freely 
placed. There are additional restrictions on their posi¬ 


tions, whose number is equal to the number of boundary 
components minus one. These restrictions will be dis¬ 
cussed in our forthcoming paper. Moreover, there is a 
“soft” constraint on the number of vortices and antivor¬ 
tices in relation to the connectivity of the system, present 
in or near the equilibrium. It is of the same “soft” nature 
as the topological solitons themselves in finite systems. 
Really, only in the infinite 2-d magnet the topological 
charge is conserved and there is an infinitely high energy 
barrier^ separating the states with different topological 
charge. In finite magnets this barrier is finite and vortices 
may enter/exit the particle through the boundary, chang¬ 
ing the topological charge of the magnetization texture 
inside. Yet, once this barrier is holding (and the magneti¬ 
zation vector is pinned to the boundary of the nanomag¬ 
net) the vortices inside the planar nanoelements behave 
like true topological solitons. 

Complete solution of the micromagnetic problems, us¬ 
ing these parametrized magnetization textures as trial 
functions, will require computation and minimization of 
the total energy of the magnet (including the exchange, 
magnetostatic and other energy terms), or, solution of 
the equations of the magnetization dynamics using the 
parameters as collective variables^ In many cases this 
can be a very complex task, which is probably still easier 
to approach numerically (as it was done in the frame¬ 
work of the Magnetism@home distributed computational 
project^). Yet, even without the total energy computa¬ 
tion, these analytical trial functions can still be useful to 
understand, classify and interpret the magnetization tex¬ 
tures, obtained in the experiment and simulations. Also, 
they can lead to beautiful analytical results, with the 
generality well beyond what numerical micromagnetics 
may offer. Like, for example, the formula for the mag¬ 
netic vortex precession frequency^ based on the dis¬ 
placed vortex model ([3]). The presented “complex vari¬ 
ables” approach might generate more useful models, now 
in multiply-connected case as well. 
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